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AB - In this filter, optimum filter coefficients are determined from 
widely scattered measurement values, for example bearings, in an 
iterative signal processing circuit. The filter coefficients and 
the measurement values are used for driving an accumulation circuit 
which determines correction values for an initial estimate of the 
target data. According to the invention, the filter is improved bv 

additional multiplication circuits for determining the coefficients 
for the dynamic target data, that is to say. the velocity 
components. 

Such filters are used in solving navigation tasks, particularly the 
passive determination of position and motion data of target 
vehicles . 
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@ Filter zur Bestimmung von Zieldaten 

Es wird ein Fitter zur Bestimmung von Zieldaten angege- . 
ben, bei dem aus stark streuenden Me&werten, z. B. Peilun- 
gen, in einer iterativen Signalverarbeitungsschaitung opti- 
male Filterkoeffizienten ermittelt werden. Mit den Filterkoef- 
fizienten und den MeBwerton wird eine Akkumuf ationsschal- 
tung angesteuert, die Korrekturwerte fur eine Anfangs- 
schatzung der Zieldaten bestimmt. Das Filter wird erfin- 
dungsgemaE durch zusatzliche Multiplikationsschaitungen 
zum Ermrtteln der Koeffizienten fur die dynamischen Zielda- 
ten, d. h. die Geschwindigkeitskomponenten, verbessert. 
Derartige Filter werden bei der Losung von Navigationsauf- 
gaben, tnsbesondere der passiven Bestimmung von Posi- 
tions- und Bewegungsdaten von Zielfahrzeugen angewen- 
det. 
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Patehtanspruche 

1. Filter zur Bestimmung von Zieldaten aus aufgenommenen Meflwerten (BBm& z. B. Peilungen, mit eraera 
Zielvektorspeicher (10) fQr einen Eingangsvektor #Jt der Zielpositionskoordinaten (fa Vektorkom- 

5 ponenten aufweist, mit einem Eingabespeicher (60), der Koordinaten jeweils eingenommener (X& Ye) 

enthalt, mit einem mit dem Zielvektorspeicher (10) und dem Eingabespeicher (60) verbundenen Koeffizien- 
tenrechner (50) zum Bestimmen einer jedem MeBwert (B^) zugeordneten Folge von MeBkoeffizienten (h$\ % 
hi?) und einer MeBwertschatzung aus den Eingangsvektor (P) und der zugehdrigen MeBpositionen 
(X& YeO> mit einer den MeBkoeffizientenausgangen des Koeffizientenrechners (50) nachgeschalteten Re- 

io chenvorrichtung (80) f ttr Filterkoeffizienten (f, u f i7 )> in der die Matrix (F) der Filterkoeffizienten (fit, in) aus 

der Matrix (H)der MeBkoeffizienten (A, u A/?)gemaB der Matrixgleichung 

F — (H' x H)~ l • H' 

15 bestimmt wird, mit einer zum Bestimmen eines Ausgangsvektor (Ap)dcr Zieldaten ausgebildeten Akkumu- 

lationsschaltung (88), in der f Or jede Vektorkomponente des Ausgangsvektors (Ap) die mit den Filterkoeffi- 
zienten (fn. fa) gewichteten Differenzen (Smi) zwischen dem MeBwert (BmO und der jeweiligen MeBwert- 
schatzung (Bhfi) komponentenweise aufsummiert werden, und mit einer Vergleicherschaltung (90) fQr den 
Ausgangsvektor {Ap), die einerseits ein Einschreiben des um den Ausgangsvektor (Ap) korrigierten Ein- 

20 gangsvektors (p) in den Zielvektorspeicher (10) bewirkt und die andererseits bei tJhterschreiten einer 

vorgebbaren Schwelle eine Ausgabe des korrigierten Eingangsverkehrs (p) oder bei Oberschreiten der 
Schwelle eine Wiederholung der Bestimmung der Zieldaten ausldst, dadurch gekennzeichnet, daB der im 
Zielvektorspeicher (10) abgespeicherte Eingangsvektor (p) ais weitere Vektorkomponente Geschwindig- 
keitskomponenten (9» Vy> V, i, V, j) aufweist, daB ein mit einem MeBzeitgeber, der Teil einer Steuervorrich- 

25 tung (77) ist, verbundener MeBzeitenspeicher (70) vorgesehen ist, daB fOr jede Geschwindigkeitskomponen- 

te (V* Vj> Vxu VrO ein Multiplizierer (25.1, 25.2, 253, 253) vorgesehen ist, der eingangsseitig mit dem 
Zielvektorspeicher (10) und dem MeBzeitenspeicher (70) verbunden ist, daB einem Zielpositionskoordina- 
ten-Eingang (51.1, 52.1) des Koeffizientenrechners (50) jeweils ein Addierer(26.1, 26.2) vorgeschaltet ist, der 
eingangsseitig mit dem Multiplizierer (25.1, 25.2, 253, 25.4) und mit dem Zielvektorspeicher (10) zur 

30 Obernahme der Zielpositionskoordinaten ' (fa fo) verbunden ist, daB an jedem MeBkoeffizientenausgang 

(56.1, 56.2) des Koeffizientenrechners (50) jeweils mindestens eine mit dem MeBzeitenspeicher (70) verbun- 
dene Multiplizierstufe (393, 39.4, 393, 39.6) zum Bilden weiterer MeBkoeffizienten (A/3, A/4, A/* A, 6 ) ange- 
schlossen ist und daB die Multiplizierstufen (393, 39.4, 393, 39.6) ausgangsseitig mit der Rechenvorrichtung 
(80) verbunden sind 

35 2. Rlter nach Anspruch 1, dadurch gekennzeichnet, daB der Koef fizientenrechner (50) eine Abstandsrechen- 

schaltung (40) aufweist, die mit Zielpositionskoordinaten (fa fa) und MeBpositionskoorcBnaten (X& Ye) 
ansteuerbar ist und Ausg&nge (46.1, 46,2) fttr die Abstandskomponenten (ky* und einen Abstandsqua- 
dratausgang (47) fQr die Summe der Quadrate der Abstandskomponenten aufweist und daB der 

Koeffizientenrechner (50) derart ausgebildet ist, daB die Bestimmung der MeBkoeffizienten (ha, ha) gem&B 

40 

hn ~ * Ha = 

45 und die Bestimmung der MeBwertschatzung (& M j) als Funktion der Abstandskomponenten, vorzugsweise 

gem&fi 

4*-arctg4*, 

50 R * 

erfolgt 

3. Filter nach Anspruch 1 oder 2, dadurch gekennzeichnet, daB die Rechenvorrichtung (80) zusatzlich mit 
Stfltzkoeffizienten (h s \, h f2t A, 3, A, 4 ) einer StOtzkoeffizientenschaltung (120) ansteuerbar ist, wobei die 

55 Matrix (F) der Filterkoeffizienten (in bis Z/4. ist bis & 4 ) aus der um die Stfltzkoeffizienten (hsu hs* hs* hst) 

erweiterten Matrix (H)der MeBkoeffizienten (h iu A/2, A, 3, A, 4 ) bestimmt wird, und daB der Akkumulations- 
schaltung (88) die Differenz (8 M $ aus den MeBwerten (Bw) und den MeBwertschatzungen BmO *>wie 
mindestens ein Dif f erenzwert (8s) aus einer Stfltzwerteingabe (S) in einen Stfitzwertspeicher (67) und einer 
von der StOtzkoeffizientenschaltung (120) bestimmten StQtzwertschatzung (S)m der gieichen Reihenfolge, 

60 in der die Rechenvorrichtung (80) mit dem MeBkoeffizienten (hn, A/ 2, A/3, A/4) und den Stfltzkoeffizienten 

(hs t, As* A53, A54) angesteuert worden ist, zugefflhn werden. 

4. Filter nach einem der Ansprflche 1 bis 3, dadurch gekennzeichnet, daB der Eingabespeicher (60) zum 
Speichem mindestens eines Stutzwertes (S), z. B. eines beobachteten Entfernungswertes (SrX eines Ge- 
schwindigkeitswertes (Sv) und/oder eines Kurswertes (SkX einer Bezugszeit (Tr) und von Beobach tungsko- 

65 ordinaten (X& Yr) ausgebildet ist und daB die StOtzkoeffizientenschaltung (120) mit ihren Eingangen (121.1, 

121.2, 122.1, X2UZ) ffir Eingangsvektorkomponenten (fa fa fa V y ) mit dem Zielvektorspeicher (10) und mit 
ihren Eingangen (123, 1213, 1223) fur die Bezugszeit (Tr) und die Beobachtungskoordinaten (Xr, Yr) mit 
dem Eingabespeicher (60) verbunden ist 
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4 5. Filter nach einem der Anspruche 3 bis 4, dadurch gekennzeichent, daB die StQtzkoeffizientenschaltung 
* (120) einen StQtzkoeffizientenrechner (150) aufweist, dessen Eing&nge (151.1, 152.1) fttr Zielpositionskoordi- 

naten in gleicher Weise wie der Koeffizientenrechner (50) Qber Multiplizerer (125.1, 125.2) und Addierer 
(126.1, 126.2) ansteuerbar sind, wobei der MultipHzierer (125.1, 125.2) Qber den Beobachtungszeiteingang 
(123) mit dem Eingabespeicher (60) verbunden ist, und dessen Eing&nge (151 .2, 152.2) fQr Mefipositionskoor- 5 
dinaten mit den Teilspeichern (63, 64) fOr Beobachtungspositionen des Eingabespeichers (60) verbunden 
sind, daB der StQtzkoeffizientenrechner (150) AusgSnge (156.1, 156.2) fur erste StQtzkoeffizienten (hsu (hsi) 
und einen Ausgang (157) fur die StQtzwertschatzung (§) aufweist, daB die StQtzkoeffizientenschaltung (120) 
Multiplizierstufen (1393, 139.4) aufweist, die eingangsseitig mit den Ausgangen (156.1, 156^) fflr erste 
StQtzkoeffizienten (hsu A52) und mit dem Beobachtungszeiteingang (123) verbunden sind, an denen aus- 10 
gangsseitig weitere Sttitzkoeffizienten (hs$, As 4) anstehen, die zusammen mit den ersten Stutzkoeffizienten 
(hs u h s 2) A usgangssignal e der StQtzkoeffizientenschaltung (1 20) bildea 

6. Filter nach einem der AnsprQche 3 bis 5, dadurch gekennzeichnet, daB den Ausgangen (156.1, 156.2) des 
Stfltzkoeffizientenrechners (150) zum Bewerten der ersten StutzkoefHzienten (hsu hsi) und der Differenz- 
schaltung (185) jeweils Gewichtsmultiplizierer (135.1, 135.2, 183) nachgeschaltet sind und daB die Gewichts- 15 

• mukiplizierer (135.1, 135.2, 183) mit einem Gewichtsfaktorspeicher (66) des Eingabespeichers (60) verbun- 
den sind. 

7. Filter nach einem der Anspruche 3 bis 6, dadurch gekennzeichnet, daB die StQtzkoeffizientenschaltung 
(120) Multipliziereinheiten (124.1, 124.2, 1243, 124.4, 139.1, 139.2) aufweist* die einerseits Qber einen Steuer- 
eingang (128) mit dem Eingabespeicher (60) und andererseits jeweils mit den Eing&ngen (121.1, 122.1) fflr 20 
Zielpositionskoordinaten (To, fb) sowie den EingSngen (1213, 1223) fflr Beobachtungskoordinaten (X& Yr) 
bzw. den Ausgangen f Or erste StQtzkoef f izienten (hs 1, A52) nachgeschaltet sind. — 

8. Filter nach einem der AnsprQche 3 bis 7, dadurch gekennzeichnet, daB die StQtzwerteingabe (S) erne 
Entfernungsstutzung (Sr) ist, daB die Sttttzkoeffizientrechner (150) die Abstandsrechenschaltung (40) auf- 
weist, deren Abstandskomponentenausgange (46.1, 46.2) mit den Ausgangen (156.2, 156.1) des Stfltzkoeffi- 25 
zientenrechners (150) fflr erste Stfltzkoeffizienten (hsu hsi) und deren Abstandsquadratausgang (47) mit 
dem Ausgang (157) fQr die StQtzwertschatzung (^ verbunden ist, und daB im Steuerspeicher (68) eine Eins 
gespeichert ist. 

9. Filter nach einem der AnsprQche 3 bis 7, dadurch gekennzeichnet, daB die StQtzwerteingabe (S) eine 
GeschwindigkeitsstQtzung (Sv) ist, daB der StQtzkoeffizientenrechner (150) die Abstandsrechenschaltung 30 
(40) aufweist, deren Abstandskomponentenausgange (464, 46.2) mit den Ausgangen (156.2, 156.1) des 
Stfltzkoeffizientenrechners fQr erste StQtzkoeffizienten (hsu /fc2)und deren Abstandsquadratausgang (47) 
mit dem Ausgang (157) des Stfltzkoeffizientenrechners (150) fflr die StQtzwertschatzung (S) verbunden sind, 
und daB im Steuerspeicher (68) eine Null gespeichert ist 

10. Filter nach einem der AnsprQche 3 bis 7, dadurch gekennzeichnet, daB die StQtzwerteingabe (S) eine 35 
KursstQtzung (Sk) ist, daB die StQtzkoeffizientenschaltung (120) als StQtzkoeffizientenrechner (ISO) dem 
Koeffizientenrechner (50) aufweist und daB im Steuerspeicher (68) eine Null gespeichert ist 

11. Filter nach einem der Anspruche 1 bis 9, dadurch gekennzeichnet, daB die Rechenvorrichtung (80) zum 
Bestimmen der Filterkoeffizienten (f) ein Feld- oder Matrixprozessor ist 

12. Filter nach einem der AnsprQche 1 bis 1 1, dadurch gekennzeichnet, daB die Steuervorrichtung (77) zur 40 
Taktgenerierung und Synchronisation und der zugehdrige MeBzeitgeber zum Erzeugen der MeBzeiten 
(Tmi) vorgesehen sind 

13. Filter nach einem der AnsprQche 1 bis 12, dadurch gekennzeichnet, daB ein mit der Peilvorrichtung (87) 
und der Steuervorrichtung (77) verbundener Manoverdetektor (76) und ein Mandverzeitspeicher (71) 
vorgesehen sind, in dem eine Anzahl (j) von Mandverzeiten (Tao, Tax) eines in aufeinanderfolgenden 45 
Zeitintervallen unterschiediiche Geschwindigkeitskomponenten (P* V x 1, V r i) aufweisenden Ziels abge- 
speichert sind, daB zwischen dem MeBzeitenspeicher(70) einerseits sowie den Multiplizierern (25.1 des 25.4) 
und den Multiplizierstufen (393 bis 39.6) andererseits eine Zeitvergleichsschaltung (79) eingeschaltet ist, die 
eingangsseitig mit dem Mandverzeitspeicher (71) verbunden ist und zeitintervalispezifische Zeitausgange 
(79.0, 79.1) fQr aus Differenzen von MeB- und Mandverzeiten gebildete Zeitsignalwerte (Tm Tmx) auf- 50 
weist, daB die Mukiplizierer (25.1, 25.2 bzw: 253, 25.4) fur Geschwindigkeitskomponenten (V* ^bzw. V xXi 
V y \) des gleichen Zeitintervalls mit dem gleichen zeitintervallspezifischen Zeitausgang (79.0, 79.1) verbun- 
den sind, und daB die jedem MeBkoeffizientenausgang des Koeffizientenrechners (50) entsprechend der 
Anzahl der Zeitintervalle nachgeschaltete Multiplizierstufen (393 bis 39.6) jeweils mit dem zugehorigen 
Zeitintervallspezifischen Zeitausgang (79.0, 79.1) verbunden sind. 55 

14. Filter nach Anspruch 13, dadurch gekennzeichnet, daB die Zeitvergleichsschaltung (79) derart ausgebil- 
det ist, daB an dem zeitintervallspeziFischeh Zeitausgang (79.0, 79.1) als Zeitsignalwert (Fia T m Tmx) 
entweder ein Nullsignalwert ansteht, wenn die MeBzeit (T M j)k\emtr oder gleich der Manoverzeit (T A o» Ta 1) 
der unteren Intervallgrenze ist, oder ein Diff erenzsignalwert zwischen MeBzeit (T yj) und Manoverzeit (Ta 0, 

T A 1) der unteren Intervallgrenze ansteht, wenn die MeBzeit (Tm) groBer als die untere und kleiner oder 60 
gleich der Mandverzeit (Tao, Tax) der oberen Intervallgrenze ist, oder ein Diff erenzsignalwert zwischen 
oberer und unterer Intervallgrenze ansteht, wenn die MeBzeit (T^groBer als die Manoverzeit (Ta x) der 
oberen Intervallgrenze ist 

15. Filter nach einem der Anspruche 13 oder 14, dadurch gekennzeichnet, daB der Mandverzeitspeicher (71) 
Qber den Steuerbus (78) mit der Rechenvorrichtung (80) ist und daB die Rechenvorrichtung (80) derart 65 
ausgebildet ist, daB der Grad der Matrizen JF/und F proportional derjenigen Anzahl (j) von Mandverzeiten 
(Ti/Jbegrenzbar ist, die bis zur jeweils aktuellen MeBzeit (TM^detektiert worden ist 

16. Filter nach einem der AnsprQche 13 bis 15, dadurch gekennzeichnet daB der Mandverzeitspeicher (71) 

[ 3 
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einen mit dem Steuerbus (78) verbundenen Zahler (713) zum Bestimmen der Anzahl Q) detektierten 
Mandverzeiten (T^und vom Zahler (713) adressierbare Speicherpiatze (71.0, 71.1) fGr die Mandverzeiten 

l^ilta^ach einem der Anspruche 1 bis 16, dadurch gekennzeichnet daB in der Vergleichseinrichtung (90) 
die Vektornorm \\Ap\\ des Ausgangsvektors (Ap) 

i Apt « Tu^o + 5> K • n u ? + td ft + 2> ^ ■ tu?T 

L > y 

mit der vorgebbaren Schwelle verglichen wird 

Beschreibung 

15 Die Erfindung betrifft ein Filter zur Bestimmung von Zieldaten der im Oberbegriff des Anspmchs 1 angegebe- 
nen Art 

Derartige Filter sind bei der L6sung von Navigationsaufgaben erforderlich, um aus mehrfach vorliegenden 
MeBwerten die Position eines zu vermessenden Objektes, beispielsweise einer Landmarke oder eines Zielfahr- 
zeugs, zu bestimmen. Insbesondere Unterwasserf ahrzeuge benutzen zur Positionsbestimmung von Zielfahrzeu- 

20 gen die passive akustische Peilung von Gerauschquellen, um den eigenen Standort durch aktive Sendesignale 
nicht zu verraten; sie sind deshalb in hohem MaBe auf solche Filter angewiesen, mit denen aus passiven 
Peilungen Positionskoordinaten bestimmt werden kdnnen, 

Es ist ein Filter, d- h. ein einer MeBvorrichtung nachzuschaltendes Signalverarbeitungssystem, bekannt, mit 
dem aus den stark streuenden MeBwerten Ausgangswerte ermittelt werden, die mit einer bekannten oder 

25 bestimmbaren statischen Sicherheit Positionen der vermessenen Objekte sind. Durch das Filter wird dabei die 
Glattung streuender MeBwerte, z. B. Peilungen, und die Umsetzung der MeBwerte in Zustande des Objektes, 
z.B. Positionskoordinaten, geleistet Diese Zustande werden auch als die Komponenten eines Zielvektors 
bezeichnet 

Zur Realisierung dieses Filters dient ein Signalverarbeitungssystem, bei dem eine MeBwertetnsehrankung und 

30 die erforderlichen Filterkoeffizienten aus einem geschatzten Eingangsvektor, namlich den Zielpositionskoordi- 
naten, und der zugehorigen MeBposition ermittelt werden. Durch das Signalverarbeitungssystem, das einem 
Modell fflr die Umsetzung der MeBwerte in die Objektzustande entspricht. wird der Fehler zwischen dem 
MeBwert und der MeBwertschatzung minimiert, so daB als eigentlicher Filter-Eingangswert die Differenz 
zwischen MeBwert und MeBwertschatzung mit dem jeweiligen Filterkoeffizienten multipliziert wird Die Sum- 

35 me aller Multiplikationsergebnisse ergibt den Wert der jeweiligen Komponente des Ausgangsvektors, in diesem 
Fall einen Fehlervektor der Positionskoordinaten, um den Eingangsvektor zu korrigieren. Der so bestimmte 
Ausgangsvektor wird komponentenweise einer Vergleichseinrichtung zugefuhrt. Weicht der Ausgangsvektor 
von einer vorgegebenen Schwelle ab, so wird der um den Ausgangsvektors korrigierte Eingangsvektor dem 
Signalverarbeitungssystem in der Art einer Iterationsschleife erneut zugeftthrt Die mehrfach iterative Filterung 

40 der MeBwerte ergibt dann letztendlich die Zielzustande, also hier die Zielpositionskoordinaten, mit einer 
Fehlerabweichung, die komponentenweise kleiner als die vorgegebene Schwelle ist 

Ein solcher Filter ist jedoch fQr die Verarbeitung von Zielzustanden ungeeignet, die sic* fOr bewegte Zide 
ergeben. Die weiteren Zieldaten , Kurs und Geschwindigkeit, werden in keiner Form berQcksichtigt Das Filter 
versagt sogar bereits nach wenigen Iterationen auch bei der Bestimmung der Zielpositionskoordinaten, da der 

45 Fehlervektor nicht mehr konvergiert 

Der Erfindung liegt die Aufgabe zugrunde, ein Filter der eingangs genannten Art so zu verbessern, daB es mit 
geringem Schaltungsaufwand die Zieldaten auch von bewegten Zielen mit mdglichst wenigen Iterationen 
bestimmt 

Diese Aufgabe wird bei einem Filter der im Oberbegriff des Anspruchs 1 definierten Art erfmdugnsgemaB 

so durch die im Kennzeichenteil des Anspruchs 1 angegebenen Merkmale geldst 

Dazu ist bei dem erfindungsgemaBen Filter der Zielvektorspeicher erweitert, um audi die die Bewegung des 
Ziels nach Geschwindigkeit und Kurs beschreibenden Geschwindigkeitskomponenten aufzunehmen. Aus den 
Zielpositionskoordinaten zu Beginn der Messung, auch als statische Zielpositionskoordinaten bezeichnet den 
Geschwindigkeitskomponenten und der jeweiligen MeBzeit werden mittels Multiplizierer und Addierer sog. 

55 dynamische Zielpositionskoordinaten gebildet, mit denen der Koeffizientenrechner angesteuert wird Aus den 
vom Koeffizientenrechner ermittelten ersten MeBkoeffizienten, dem sog. statischen Teil der MeBkoeffizienten, 
werden in den nachgeschalteten Multiplizierstufen durch erneute Multipiikation mit der MeBzeit nunmehr 
weitere MeBkoeffizienten, der sog. dynamische Teil der MeBkoeffizienten, gebildet In einer Rechenvorrichtung 
fQr Matrixoperationen erfolgt die Umwandlung der MeBkoeffizienten in die erforderlichen Filterkoeffizienten, 

60 mit denen dann in komponentenweise gleichartiger Verarbeitung die Zieldaten besdmmt werden. 

Das erfindugnsgernaBe Filter hat den Vorteil, daB zum Bestimmen des dynamischen Teils der MeBkoeffizien- 
ten die gleichen Schaltungen, Rechenvorrichtungen, wie fQr die statischen MeBkoeffizienten benutzt werden. Es 
sind nur wenige zusatzliche Multiplizierer, Multiplikadonsstufen notwendig, um die Zeitabhangigkeit der dyna- 
mischen Koeffizienten zu erfassen. Fur die nachgeschaltete Rechenvorrichtung, die Akkumulationsschakung, 

65 den Vergleicher usw. werden entsprechend der Anzahl der berucksichtigten Vektorkomponenten jedoch Ver- 
vielfachungen der intemen Verarbeitungsschaltung notwendig. Bei geringen Anforderungen an die Verarbei- 
tungszeit sind daher auch im Multiplex genutzte Verarbeitungsschaltungen, d h. "Komponenten-Filter", reali- 
sierbar. 
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Weitere vorteilhaf te Ausgestaltungen der Erfindung ergeben sich aus den UnteransprQchen. 

So ist es insbesondere vorteilhaft, die Erfindung gem&B dem Anspruch 3 weiterzubilden, da durch die 
Ermittlung von zusatzlichen Stutzkoeffizienten mit einer Stfltzkoeffizientenschaltung, die im wesentiichen der 
MeBkoeffizientenschaltung zum Bestimmen der MeBkoeffizienten entspricht, dem Filter StOtzwerte eingegeben 
werden kdnnen, die durch Beobachtung oder andere MeBverfahren ermittelt worden sind. Die Stfltzkoeffizien- 5 
ten und die zugeh6rige Stfltzwertschatzung werden aus den Vektorkomponenten des Eingangsvektors und 
erforderlichenfalls zusfitzlichen StQtzwerteingaben bestimmt Die MeBkoeffizientenmatrix //wird zusatzlich mit 
den Stutzkoeffizienten erweitert und infolgedessen wird auch eine zur Ansteuerung der Akkumulatorschaltung 
erweiterte Filter-Koeffizientenmatrix F gebildet Die Folge davon ist, daB die Filterung der MeBwerte ein- 
schlieBIich der zusfctzlich in der Akkumulationsschaltung verarbeiteten Differenzwerte schneller, d h. mit einer / 10 
wesentiichen geringeren Zahl von iterierenden Systemdurchiaufen, zum genauen Zieldatenvektor filhrt 

Es ist weiterhin vorteilhaft gemaB Anspruch 6 die Differenzwerte und die Stutzwertkoeffizienten mit einem 
Gewichtsfaktor zu bewerten. Gute StiltzWerteingaben, Z.B. Radarmessungen, werden mit groBem Gewicht 
berucksichtigt, so daB ein schnelles Einschwingen des Filters auf die wahren Werte des Zieldatenvektors erfolgt 
Weniger gute StOtzwerte, z. B. ungenauere Sehrohrbeobachtungen, werden mit geringerem Gewicht berflck- 15 
sichtigt und verandern dann das Einschwingverhalten kaum. Auf diese Weise ist das Filter der jeweiligen MeB- 
und Beobachtungssituation optimal anpaBbar. 

Derartige Filter zur Bestimmung von Zieldaten durch Filterung passiver Peilungsmessungen konvergieren 
erst dann zu einem stabilen Ergebnis, wenn das messende Fahrzeug selbst ein Mandver; das sog. Eigenmandver, 
gefahren hat Hierbei ist die Beriicksichtigung von StOtzwerten zum Bestimmen der erweiterten Filterkoeffi- 20 
zientenmatrix Fbesonders vorteilhaft, da das Eigenmandver vermieden werden kann, wenn sehr gute Messun- 
gen der Stfltzwerte mit hoher Gewichtung ausgewertet werden. 

Der erfindungsgemaBe Aufbau des Filters weist daruber hinaus den Vorteil auf, daB die Beriicksichtigung von 
Stutzwerten fur jede Iteration unterschiedlich gesteuert werden kann, da in jeder Iteration die Filterkoeffizien- 
tenmatrix neu berechnet wird 25 

Die Erfindung ist anhand eines in der Zeichnung dargestellten- Ausfuhrungsbeispiels im folgenden naher 
beschrieben. Es zeigt 

Fig. 1 ein Blockschaltbild einer Vorrichtung zur Bestimmung von Zieldaten, 

Fig. 2 Ein Blockschaltbild einer Koeffizientenschaltung, 

Fig* 3 ein Blockschaltbild einer Abstandsrechenschaltung, 30 
Fig. 4 ein Blockschaltbild einer Stfltzkoeffizientenschaltung, 

Fig. 5 ein Blockschaltbild einer Vorrichtung zur Bestimmung der Zieldaten eines Mandver fahreriden Ziels, 
Fig. 6 ein Blockschaltbild einer Koeffizientenschaltung zur MeBkoeffizientenbestimmung bei einem Mandver 
f ahrenden Ziei, 

Fig. 7 eine Darstellung einer Zielbahn eines Mandver fahrenden Ziels. 35 

Fig, 1 zeigt das Blockschaltbild einer Vorrichtung mit einem Filter zum Bestimmen von Zieldaten aus 
MeBwerten B Mh Die Vektorkomponenten J$b, fb, V* % dieser Zieldaten sind als Eingangsvektor p in einem 
Zieldetektorspeicher 10 mit Speicherelementen 1 1, 12 fur Zielpositionskoordinaten X& fb und Speicherelemen- 
ten 13, 14 fQr Geschwindigkeitskomponenten Vx, t^abgelegt Dem Zieivektorspeicher 10 ist eine Koeffizienten- 
schaltung 20 nachgeschaltet, an die eingangsseitig ein Eingabespeicher 60 mit seinen Teilspeichern 61, 62 fur 40 
MeBpositionen X& Ye und ein MeBzeitenspeicher 70 angeschlossen sind Die Koeffizientenschaltung 20 ist zur 
Obergabe von MeBkoeffizienten A/i, A/2, A;* A/4 mit einer Rechenvorrichtung 80 und zur Obergabe einer 
MeBwertschatzung 4m? mit einer Differenzeinheit 85 verbunden. Die Differenzeinheit 85 ist eingangsseitig auch 
an einen MeBwertspeicher 86 angeschlossen, in dem die von einer Peilvorrichtung 87 zu den MeBzeiten Tat, 
auf genommenen Peilungen als MeBwerte Bw abgespeichert sind 45 

Die Rechenvorrichtung 80 ist derart ausgebildet, daB die vollstandige Matrix //der MeBkoeffizienten hiu A/2, 
A/3, A/4 gespeichert und in eine Matrix Fder Filterkoeffizienten fa t fa fa fa umgerechnet wird, mit denen eine 
Akkumulationsschaltung 88 ansteuerbar ist 

Die Akkumulationsschaltung 88 ist auBerdem eingangsseitig mit der Differenzeinheit 85 verbunden, so daB 
jeder Filterkoeffizient fii, fa fa fa mit der zugehorigen Differenz 6 M i aus dem MeBwert B Mi und der MeBwert- 50 
schatzung 6 M i multipliziert und aufsummiert wird und einen Ausgangsvektor Ap der Akkumumulationsschal- 
tung 88 bildet Ausgangsseitig ist die Akkumulationsschaltung 88 mit einer Vergleichsschaltung 90 und flber eine 
an dem Zieivektorspeicher 10 angeschlossene Additionsstufe 94 mit einer Torschaltung 95 verbunden, fiber die 
der Zieivektorspeicher 10 mit dem urn den Ausgangsvektor Ap der Akkumulationsschaltung 88 korrigierten 
Eingangsvektor p ansteuerbar ist Die Vergleichsschaltung 90 ist ausgangsseitig an die Torschaltung 95 und an 55 
ein Ausgabe-Interface 96 zum Ausgeben des Eingangsvektors p t & h, der Zieldaten, angeschlossen. Die Zusam- 
menschaltung des Zielvektorspeichers 10 und des Ausgabe-interface 96 erfolgt flber eine Verbindungsleitung97, 
auf die auch die Additionsstufe 94 zugreift 

Der Zieivektorspeicher 10 ist jedoch nicht nur uber die Torschaltung 95 mit dem urn den Ausgangsvektor Ap 
korrigierten Eingangsvektor p % sondern auch von einem Basisdatenspeicher 15 fQr Vektorkomponenten Xo, Yq, so 
V n (^ansteuerbar. 

Bei der Bestimmung der Zieldaten wird von Zielen ausgegangen, die sich abschnittsweise^gleichf6rmig 
bewegen. Das Ziel fahrt also unbeschleunigt, d h. mit konstanter Geschwindigkeit, auf konstantem Kurs. Diese 
konstante Geschwindigkeit und der Kurs werden durch das Filter ermittelt, d. h., die von der Geschwindigkeit 
und der Zeit abhSngigen Filterkoeffizienten fa fa fa fa werden bestimmt und die Peilungen so gefiltert, daB 65 
auBer den Zielpositionskoordinaten Xo, % auch die Geschwindigkeitskomponenten V* *> ausgegeben werden, 
aus denen sich aktuelle Schatzwerte f flr Entf ernung, Kurs und Geschwindigkeit des Ziels ergeben. 

Das Filter weist dazu einen Zieivektorspeicher 10 fiir die Zieldaten auf, in dem zusatzlich zu den statischen 

5 
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Zielpositionskoordinaten X* ?o auch die Geschwindigkeitskomponenten V» % abgelegt sind. Da bei passiven 
Ortungsvorrichtugnen fttr fremde Ziele nie die wirklichen sondern nur sog. geschatzte Zieldaten als Zielvektor 
P verfttgbar sind, sind die Komponenten durch ein zusatzliches dachfdrmiges Zeichen ( ~) kenndich gemacht 
Die dabei auf tretenden Fehler AX, A Y, A V* A V y geben dann die Unterschiede zwischen den geschatzten und 
den wirklichen Zieldaten an. b 

Das Filter ist dasjenige Obertragungssystem, mit dem aus einer Vielzahl /« 1 bis N gemessener Peilungen 3m 
der Vektop p der Zieldaten Xo, *o, V* Yy 

C ™us*den statischen Zielpositionskoordinaten Xb, ?o, den Geschwindigkeitskomponenten und alien 

gespeicherten MeBzeiten Tm zu denen die Peilungen £mi gemessen worden sind, werden nacheinander kompo- 
nentenweise durch Multiplizierer 25.1, 25.2 und Addierer 26.1. 2fc2 dynamische Zielpositionskoordinaten 
£i- Xo+V x ■ 7m? und ?i - ?i+ % • 7m/ erzeugt und einem Koeffizientenrechner 50 zugeftthrt Der Xoeffi- 
zientrechner 50 bestimmt auf der Grundlage der MeBgeometrie fttr jeweils eine Zielposition X h R unter 
Berflcksichtigung der dazugehdrigen, gespeicherten MeBposition X&> rs eine sog. MeBwertschatzung Bm> A h. 
die Schfitzung einer Peilung von der MeBposition X& Ye zur geschatzten dynamischen Zielposition X* Y- t 

Der Koeffizientenrechner 50 bildet ein Filtersystem nach, bei dem aus einer einzigen statischen Zielposition 
X, fund der zugeh6rigen MeBposition X& Yedle geschatzten Peilungen /found erste MeBkoeffizienten hi und 
hi bestimmt werden. . 

Wird das auf die vorgenannte Weise fttr statische Zustande realisierte MeBmodell auf die dynamischen 
Zielpositionskoordinaten X b ft angewendet, so ergibt sich die aus den geschatzten Zielzustanden einschlieBhch 
der Geschwindigkeitskomponenten ermittelte, geschatzte Peilung Bmzu 

B Mi - arctg ■ : , 

oder mit den Abkttrzungen f Qr die Abstandskomponenten Axi und £ w 



jS Af/ = arctg 



mit 

wobei /P/das Abstandsquadrat zwischen der Zielposition fiund der MeBposition X& ft und der MeBposition 
Xek ybdarstellt . - . , v# ot 

Mit den vom Koeffizientenrechner 50 aufgrund der dynamischen Zielposition X b y/ermittelten ersten MaBko- 
effizienten hn und hn stehen bereits die ersten Elemente einer MeBkoeffizientenfolge fest, deren Elementenan- 
zahl von der Komponentenanzahl des Zielverstarkers p abhangt Weitere zur Bestimmung des vollstandigen 
Zielvektors p notwendige MeBkoeffizienten A l3 und h lA werden zusatzlich durch Multiplikation der ersten 
MeBkoeffizienten hn und Zr/ 2 mit der MeBzeit 7A#?ermittelt 
Zum Nachweis dieser Zusammenhange wird in einer der GL (1) entsprechenden Beziehung fur die gemessene 
so Peilung Bm der wirkliche Zielzustandsvektor durch den geschatzten Zielzustandsvektor p und einen Fehiervek- 
tor Ap 

Ap - (AX^AYo^AV^AVy)' ( 4 ) 

ersetzt Aus der Taylor-Reihenentwiklung dieser Beziehung nach dem Fehlervektor Ap ergibt sich durch 
Weglassen der Glieder hoherer Ordnung ein linearisiertes MeBmodell fttr die dynamische Zielzustande zu 

B Mi -B Mi -^AX,^^AY 0 -^ T Mr AV x ^^ T Mr AV y (5) 

. mit den MeBkoeffizienten 



65 h lx = - j* , h l2 = , h i} = h f \ ■ T MI und h i4 = h n ■ T Ui . 
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FaBt man die MeBkoeffizienten h;x h i4 in einer Matrix //und die Differenzen zwischen MeBwert B Mi und 

MeBwertschatzung B M i in dem Vektor Szusammen, so gibt die MeBgleichung 

8-HAp < 6 ) 

die Abhangigkeit der Peilungsdifferenz 8 M , = Bm~ B Mi vom Fehlervektor 4p und der MeBkoeffizientenmatrix 
//an. Das zu realisierende Signalverarbeitungssystem ergibt fiir den Fehlervektor Ap dann einen Wert minima- 
ler Varianz, wenn aus dem Differenzenyektor 5 und einer Filterkoeffizientenmatrix Fdurch das Signalverarbei- 
tungssystem als Ausgangsvektor der Fehlervektor Ap zu 

Ap~F*8 < 7 ) 
bestimmt wir8 und wenn fur die Filterkoeffizientenmatrix F 

F — (H' • H)- x • H' ( 8 ) 15 

gilt, wobei //'die Transponierte der MeBkoeffizientenmatrix Hist 

Die Gesamtheit aller Messungen kann somit durch das Matrix-Gleichungssystem, der sog. MeBgleichung, 
gemaB Gl. (6) ausgewertet werden. Das zu realisierende Signalverarbeitungssystem ergibt fiir den Ausgangsvek- 
tor Ap dann eine optimale Schfttzung, wenn er gemaB GL (7) und GL (8) aus dem Vektor 5 der Differenzen von 20 
MeBwerten Ao und MeBwertschatzungen Bm bestimmt wird, wobei die Matrix Fdie FilterkQeffizientenmatrix 
darstellt 

Diese Filterkoeffizientenmatrix FgemaB GL (8) wird in der Rechenvorrichtung 80 bestimmt, der die Gesamt- 
heit der MeBkoeffizientenfolgen, die MeBkoeffizientenmatrix H, zugef Ohrt wird. 

Die Rechenvorrichtung 80 ist dann mit der Akkumulationsschaltung 88 zusammengeschaltet in der kompo- 25 
nentenweise die Dif ferenzen zwischen dem MeBwert B^und der MeBwertschatzung Bms mit der Filterkoef- 
fizientenmatrix Fmultipliziert und jeweils Ober alle Produkte aufsummiert werden, so daB an der Akkumula- 
tionsschaltung 88 ausgangsseitig der Ausgangsvektor Ap f flr die Zieldaten ansteht 

Der Ausgangsvektor Ap wird komponentenweise einer mit dem MeBzeitenspeicher 70.verbundenen Ver- 
gleichsschaltung 90 zugefOhrt, in der aus den statistischen und dynamischen Komponenten eine Vektornorm des 30 
Ausgangsvektors Ap gebildet und mit einer Schwelle verglichen wird Ergibt der Vergleich eine unzul&ssige 
Abweichung, so wird ein urn den Ausgangsvektor Ap korrigierte Eingangsvektor p in den Zielvektorspeicher 10 
ubertragen, und durch das Signalverarbeitungssystem werden mit den gespeicherten MeBwerten Bm und Beob- 
achtungszeiten Tms in Iterationen verbesserte Ausgangsvektoren Ap bestimmt, bis die Vergleichsbedingungen 
erf Clllt sind. . . 35 

Der in der Vorrichtung gemaB Fig. 1 vorgesehene Eingabespeicher 60 weist ferner Teilspeicher 63, 64 fOr 
Beobachtungskoordinaten Xr, Yr einer zusatzlichen MeBposition, einen Bezugszeitspeicher 65, einen Gewichts- 
faktorspeicher 66, einen StOtzwertspeicher 67 und einen Steuerspeicher 68 auf. Zur Ansteuerung der einzelnen 
Speichereinheiten ist der Eingabespeicher 60 mit einem Eingabe-Interface 69 verbunden, an das eingangsseitig 
eine PositionsmeBeinrichtung 72, eine GeschwindigkeitsmeBeinrichtung 73, eine KompaBeinheit 74 und ein 40 
Datensichtgerat75angeschlossensind -?? 

Die Vorrichtung zum Bestimmen von Zieldaten wird durch eine Stfltzkoeffizientenschaltung 120 zum Ermit- 
teln von StOtzkoeffizienten hsu h S s> h S A vervollstandigt, die eingangsseitig mit den Speicherelementen 11 bis 14 
des Zielvektorspeichers 10, mit den Teilspeichern 63, 64 fttr Beobachtungskoordinaten Xr, Yr mit dem Zeitspei- 
cher 65, dem Gewichtsfaktorspeicher 66 und dem Steuerspeicher 68 verbunden ist Ausgangsseitig ist die 45 
StOtzkoeffizientenschaltung 120 zur Obergabe der StOtzkoeffizienten h S u hs* h S A mit der Rechenvorrichtung 80 
zusammengeschaltet Die StOtzkoeffizientenschaltung 120 ermittelt darOber hinaus eine Stutzwertschatzung S, 
die in einer der StOtzkoeffizientenschaltung 120 und dem Stutzwertspeicher 67 nachgeschalteten Differenzschal- 
tung 185 in Verbindung mit der Stutzwerteingabe S einen Differenzwert 8s ergibt Die Differenzschaltung 185 
imd die Differenzeinheit 85 sind Ober einen vom Steuerbus 78 ansteuerbaren Schalter 181 wechselweise mit der 50 
Akkumulationsschaltung 88 verbindbar, so daB entweder die Differenzen 8m aus MeBwerten B M i und MeBwert- 
schatzungen B Mi oder die Differenzwerte foaus StOtzwertschatzung Sund StQtzwerteingabe Sin der Akkumula- 
tionsschaltung 88 mittels der Filterkoeffizientenmatrix F gefiltert werden und so der Ausgangsvektor Ap 

bestimmt wird. . , _. , 

In einer zentralen Steuervorrichtung 77 werden die erforderlichen Taktsignale insbesondere zum Em- und 55 
Auslesen der Speicher 10, 60, 70u der Koeffizientenschaltungen 20, 120, der Rechenvorrichtung 80 und der 
Interfaces 69, 96 generiert und Qber den Steuerbus 78 flbertragen. 

Fig. 2 zeigt im Detail die Koeffizientenschaltung 20 mit Eingangen 21.1, 22.1 fOr Zielpositionskoordinaten Ao, 
iPb, Eingangen 21.2, 22.2 fOr Geschwindigkeitskomponenten V* % Eingangen 213, 223 fflr MeBpositionskoor- 
dinaten X&, Y E und einen MeBzeiteneingang 23 sowie mit MeBkoeffizientenausgangen 36.1 bis 35.4 fOr die eo 
MeBkoeffizienten hn bis hi a und mit einem Ausgang 37 fQr die MeBwertschatzung B Mi . Die Koeffizientenschal- 
tung 20 weist Multiplizierer 25.1 und 25.2 auf, die jeweils eingangsseitig mit dem Eingang 21.2 bzw. 22^ und 
beide mit dem MeBzeiteneingang 23 verbunden sind und denen ausgangsseitig jeweils ein Addierer 26.1 bzw. 
26.2 nachgeschaltet ist, der eingangsseitig auBerdem an dem Eingang 21.1 bzw. 22.1 angeschlossen ist Den 
Addierern 26.1, 26.2 ist ausgangsseitig ein Koeffizientenrechner 50 mit Zielpositionskoordinaten-Emgangen 65 
51.1, 52.1 nachgeschaltet, der eingangsseitig ferner mit den Eingangen 213 bzw. 223 der Koeffizientenschaltung 
20 verbunden ist Der Koeffizientenrechner 50 ist ausgangsseitig mit dem Ausgang 37 fur die MeBwertschatzung 
B Mi sowie einerseits unmittelbar mit den MeBkoeffizientenausgangen 36.1 bzw. 36.2 und andererseits Qber 
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sen ist AuBerdem sind in dem Koeffizientenrechner 50 Dividierschaltungen 54.1 bis M ^^ g ^™ ™ e ?! 
Abs^dl?echen"haltun^ 40 mit ihren AusgSngen ffir Abstandskomponenten 46.1 und 462 uui einem Ab- 
SSSSSSS^ in der Weise nachgeschaltet sind. daB an den Divid.er*haltungen 54. I^Sf 

*«•♦ Hpm Aiwffanj? 157 der Stfltzwertschatzung £ verbunden ist Den Ausgangen 4&.1 una rur aie 

«n 153.1 und 1533 angesteuert werden. Die Ausgange der Summationsschaltungen 153.1, 1533 bilden die 
AuWange 156.1. 1563 fur erstere Stfitzkoeffizienten des StfltzkoefFizientenrechners 150. rmtPfllen 

AuErirtehn der Datenverarbeitungsschaltungen und ist mit bekannten iogischen Schaltungen »eicht zu 
real^S Dauber hinaus enthah die zentrale Steuervorrichtung 77 einen MeBzeitgeber, urn MeBzeiten 7m an 

etosotehe ^7ran"f g^dTe zur automatischen Zielverfolgung detektierter Zeile ein Jf« cht ?V^^ 
^^ 87 die iniedem Zeitaugenblick Ober MeBwerte des verfolgten Zieies verfOgt, wird durch geeignete 
Ab^SSmpuke der S uTrvorSung 77 abgefragt und die als MeBwerte B M ermitteiten Peilungen werden in 
d» ^ MeBweSpeicher ^ Gleichzeitig wird die in der Steuervomchtung " P"?""*^ 

« MeBzeit Sin dem MeBzeitenspeicher 70 gespeichert. Mit der Abfrage der Peilvornchtung 87 wird auch die 
AM^g ^« ^osSn^Beinrichtung 72 Ober das Eingabe-Interface 69 ausgelast und 

Xe^ft^i^nskoordinaten Xb, Yb in die Teilspeicher 61 und 62 des Emgabespeichers 60 Obertragen. In dieser 
* W^w?r^ 

rige, zusatzlich durch den Index /gekennzeichnete Werte die Peilungen Ba* die MeBzeiten 7m, und MeBposi 

Obwohl in Fig 1 nicht naher dargestellt. konnte eine Modifizierung des Basisspeichers 15 Ober das Eingabe- 
es Interfac -69 und das DaUnsichtgerlt 75 ebenso einfach realisiert werden, mit der dann ^ST^SSSt 
sfchtgerat 75 her Anfangswerte in den Basisspeicher 15 emgeschneben werden und somit eine Filt^rung der 
Swerte B M zum Besiimmen der Zieldaten *b, ft. % durchgefOhrt werden, die von veranderbaren ggf. 
genaueren Anfangswerten ausgeht 
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Die Qbernahme von Anfangswerten erfolgt jeweils nur zu Beginn eines MeBzyklus, n&mlich dann t wenn im 
'* Zielvektorspeicher 10 keine Werte der Komponenten Xq, fo, V» V y vorliegen, die genauer als die Anfangs werte 

des Basisspeichers 15 sind. 

Der im Zielvektorspeicher 10 abgelegte Eingangsvektor p = (X6, rb, v^)' enthah zus&tzlich zu den 
statischen Zielpositionskoordinaten Ab, rb die Geschwindigkeitskomponenten Vy, aus denen unter Berttck- 5 
sichtigung der MeBzeiten rjwund der Me&positionen X& Yb die Koeffizientenschaltung 20 die MeBkoeffizien- 
ten hi \ Jij2> A/3, A/4 ermittelt Dazu sind in der Koeffizientenschaltung 20 die Multiplizierer 25.1 und 25.2 sowie die 
Addierer 26.1 und 26.2 vorgesehen, auf denen aus dem Eingangsvektor p und den MeBzeiten Tmi dynamische 
Zielpositionskoordinaten gebildet werden. In der Abstandsrechenschaltung 40 des Koeffizientenrechners 50 
werden dann aus den dynamischen Zielpositionskoordinaten unter BerGcksichtigung der MeBpositionskoordi- 10 
naten Xm und die Abstandskomponenten fix; und £w gem£B GL (3.1) und (3.2) und das Abstandsquadrat U 2 / 
gem&B Gl. (33) ermittelt Aus den Ausgangssignaien der Abstandsrechenschaltung 40 werden durch die Dividier- 
schaltungen 54.1 bis 543 und die Arcus Tangens-Rechner 55 die ersten MeBkoeffizienten Jnu A/2 und die 
Me&wertschatzung 2§m/ bestimmt Mit den Multiplikationsstufen 393 und 39.4 werden dann durch erneute 
Multipiikation mit der MeBzeit Tmi aus den ersten MeBkoeffizienten An, A/2 die weiteren MeBkoeffizienten A/3 15 
und A/4 ermittelt An der Koeffizientenschaltung 20 steht somit ausgangsseitig eine Folge von vier MeBkoeffi- 
zienten A/1, A/2, A/3, A/4 an, die Qber den Index idem entsprechenden MeBwert £a*# zugeordnet ist 

Bei der Komplexitat des Filters konnen die Filterkoeffizienten fa, i\% fa fa nicht direkt und im wesentlichen 
unmittelbar anhand desjenigen Modells gem&B GL (7) bestimmt werden, das den Eingangsvektor p bzw. den 
Fehlervektor zip als Funktion von Filterkoeffizientenmatrix Fund MeBwerten Em berechnet sondern es ist 20 
erforderlich, zunSchst die MeBkoeffizienten A/ 1, A/2, A/3, A/4 des MeBmodells zu bestimmen und dann die vollst&n- 
dige MeBkoeffizientenmatrix //in die Filterkoeffizientenmatrix Fumzurechnen. Dazu ist die Rechenvorrichtung 
80 vorgesehen, in der die nacheinander bestimmten Folgen der MeBkoeffizienten An. A/ 2 , A/a, A/4 gespeichert 
werden. Durch die Steuervorrichtung 77 wird dann die Umrechnung der MeBkoeffizientenmatrix H in die" 
Filterkoeffizientenmatrix FgemSB der Matrixgleichung (8) freigegeben, wenn die durch den Index / bestimmte. 25 
Anzahl Folgen von MeBkoeffizienten A/t bis A/ 4 gespeichert ist Derartige Umrechnungen haben in der Regel 
einen grdBeren Zeitbedarf, so daB hier eine asynchrone Verarbeitbarkeit erfolgt mit einem Start durch- die 
Steuervorrichtung 77 und einer Fertigmeldung an die Steuervorrichtung 77 durch die Rechenvorrichtung 80. Mit 
der Fertigmeldung sind dann in der Rechenvorrichtung 80 alle den MeBwerten Bmi zugeordneten Filterkoeffi- 
zienten ft 1, fa, fa fa gespeichert 30 

Zeitgleich mit der Obergabe der MeBkoeffizienten A/i bis A/4 an die Rechenvorrichtung 80 ist die in der 
Koeffizientenschaltung 20 ermittelte MeBwertschStzung Bm an die Differenzeinheit 85 Obertragen worden, in 
der die jeweiligen Differenzen dM/gebildet und gespeichert werden. 

Das eigentliche Ausgangssignal des Filters, der Ausgangs- oder Fehlervektor Ap der Zieldaten, wird kompo- 
nentenweise durch die Multipiikation der Filterkoeffizienten fa bis fa jeweils mit der Differenz 8mj und der 35 
Summation all dieser Produkte uber alle / in der Akkumulationsschaltung 88 gebildet Dieser Ausgangsvektor 
Ap wird dann in der Vergleichsschaltung 90 in der Weise geprQft, daB seine Vektornorm || Ap || 

|| Ap\\ - [(AXq + AV X • T M if + (AY 0 + AV X . T Ml )*J<* (9) 

40 

bestimmt und mit einer in der Vergleichsschaltung 90 abgelegten, fest vorgegebenen Schwelle verglichen wird 
Wird die Schwelle von der Vektornorm || Ap || gemaB GL (9) unterschritten, so wird der Ausgangsvektor Ap in 
der Additionsstuf e 94 zum Eingangsvektor p addiert, das Additionsergebnis als verbesserter Eingangsvektor p in 
den Zielvektorspeicher 10 eingeschrieben und seine Ausgabe Ober das Ausgabe-Interface 96 z. B. an das 
Datensichtgerat 75 freigeschaltet Andemfalls erfolgt lediglich das Einschreiben des korrigierten Eingangsvek- 45 
tor p in den Zielvektorspeicher 10 und ein Neustart der Bestimmung der Filterkoeffizientenmatrix Fmit den 
gleichen MeBwerten Bm> jedoch diesem bereits verbesserten Eingangsvektor p. 

Die iterative Bestimmung der Zieldaten konvergiert schneller, d h. mit weniger Iterationen, wenn in einer 
zusatzlichen StUtzkoeffizientenschaltung 120 die Stutzkoeffizienten A/t, A/2, A/3, A, 4 und eine Stutzwertschatzung 
S bestimmt werden, mit der zusatziiche Messungen oder Beobachtungen, beispielsweise Entfernungsmessung . 50 
mit Radar, Geschwindigkeitsbeobachtung aufgrund des Propellergerausches oder Kursbeobachtung 1 durch 
Sehrohrbeobachtung, durch die Rechenvorrichtung 80 in die Berechnung der Filterkoeffizientenmatrix Fmit 
einbezogen werden. 

In der der Koeffizientenschaltung 20 weitgeheiid gleichartigen Stiltzkoeffizientenschaltung 120 werden die 
einem Stutzwert S zugeordneten StGtzkoeffizienten A51 bis A54 und das Quadrat der Stutzwertschatzung S 2 55 
ermittelt und vorzugsweise im AnschluB an alle MeBkoeffizienten A/i bis A l4 der Rechenvorrichtung 80 zuge- 
f Qhrt Die Rechenvorrichtung 80 berechnet eine erweiterte Filterkoeffizientenmatrix F t und im AnschluB an die 
Akkumulation der Produkte aus Differenzen 8w und Filterkoeffizienten fa bis §4 wird fiber den Schalter 181 die 
Differenz &aus dem Quadrat des Stutzwerts S 2 und dem Quadrat der Stutzwertschatzung S 1 an die Akkumula- 
tionsschaltung 88 durchgeschaltet, mit den Filterkoeffizienten fs\ bis f$4 multipiiziert und somit ein Ausgangs- 60 
vektor Ap mit den Komponenten 



65 

\ 
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A X = £/„ • 6 m +fsi *s S/« * **< +/« ' % (10) 
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^SibS gibt /- Ndie Anzahl MeBwerte an. die bis zu einem Zeitpunkt, L a. der aktuellen MeBzeit, ausgewertet 

W Ke to F«. 4 dargestellte StQtzkoefFizientenschaltung 120 ist mit geringem Aufwand fOr verschiedene Stfltz- 
werteingaben Sbzw. deren Quadrate S 2 modifizierbar, dazu sind die Multiplikationsschaltungen 124.1 bis 124.4 
vorgesehen, die ttber den Steuereingang 128 vom Steuerspeicher 68 ansteuerbar sind. n ^ ka ,* 

Mit der Eingabe einer quadrierten Entfernungsstfltznng S 2 * und der zur Bezugszeit Tjt gehdngen Beobach- 
tungsposition Xr, Yr wird in den Steuerspeicher 68 eine Eins eingegeben und damrt das Durchschalten der 
Zielpositionskoordinaten X* fb an die Addierer 26.1. 26a und der Beobachtungsposmon X* und Y*> an den 
StOtekoeffizientenrechner 150 bewirkt sowie die Aktivierung der Ausgange 136.1 und 136.2 durch die Mult.pl - 
zierstufen 139.1 und 139.2 gewahrleistet Indem Bezugsleitspeicher65 ist die Bezugszeit T* der Entf ernungsstut- 
zung Sr gespeichert, so daB dem StOtzkoeffizierungsrechner 150 das Ausgangssignal der Adddierer 126.1 und 
126^ als beobachtete dynaroische Zielpositionskoordinaten Xo + tfmi und % + V,Tr und die Beobachtungsposi- 
tionenAjb y«zugefuhrt werden. ■ 

In dem in Fig. 4 dargestellten StOtzkoeffizientenrechner 150 ergeben sich daraus erste StOtzkoeffizienten H Sr h 

hsR>7M 

hs R , - -2&i As* - -2ft, (1 U) 



mit 

30 & - Xo+ % ' Tr-XrvbAA,- ?o+ % • T*r~ Yr, 

da die Muhiplizierstufen 139.1 und 139.2 lediglich noch eine Multiplikation mit Eins bewirken, und das Quadrat 
der StOtzwertschatzung &r zu 

35 &R-&. 
wobei 

ist Die weiteren StOtzkoeffizienten 

As* - hs# • r^und/is* = • Tr 0 1 - 2 ) 

ergeben sich als Aiisgangssignale der Muhiplizierstufen 1393, 139.4. 

Mit der Eingabe einer quadrierten GeschwindigkeitsstGtzung Sv ttber das Datensichtgerit 75 werden der 
Steuerspeicher 68 auf Null und der Bezugszeitspeicher 65 gesetzt Damit sind die Einginge 121.1, 1213 und 
122.1, 122.3 inaktiviert und die Ausgange 136.1 und 136.2 zu Null gesetzt, so daB lediglich die zur Geschwindig- 
keitsstQtzung S^gehdrige Koeffizientenfolge zu 

h Sv t-0,hsv* - 0,hsv> - + * ftundAsy* - +2 % ( 12 ) 



Fur die Eingabe einer KursstOtzung S K ware gemaB einem nicht weiter dargestellten AusfOhrongsbeispiel der 
StOtzkoeffizientenrechner 150 in der StQtzkoeffiziemenschaltung 120 durch den Koeffizientenrechner 50 iden- 
tisch zu ersetzen. Wird dann bei Eingabe der KursstOtzung S K durch das Datensichtgerat 75 der Steuerspeicher 
68 zu Null und der Bezugszeitspeicher 65 zu Eins gesetzt, so werden in der Sttttzkoeffizientenschaltung 120 die 
StOtzkoeffizienten 



60 -fr p /13V 

hs K \ = 0, h Sja - 0. ^3 - ™<* - pf v 

gebildet, die Qber die Ausgange 136.1 bis 136.4 der Rechenvorrichtung 80 der MeBkoeffizientenmatrix H 

65 zugefuhrt werden. _ , _ . ~ , . 

Bei der Filterung von MeBwerten B w unter Einbeziehung einer Entf ernungsstutzung Sr bzw. ener Gerchwin- 
digkeitsstQtzung Si^werden somit die Differenzen Ssr bzw. &vaus den quadrierten Stutzwerten S/P bzw. und 
Schatzungen bzw. $v* t d h. 

10 
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Ssr - Sr * — £v bzw. S S v = Sv> — Sv* ( 14-1) 

mit den Filterkoeffizienten f s bis f s R * bzw. fs v *h\$fs v «gewichtetWohingegen bei einer KursstQtzung die 
Dif f erenz 8sk aus Unearer Kurssttl tzung Sk und -schatzung §k 

5sk-Sk-Sk (14^) 
gefiitert wird. 

Das AusfQhrungsbeispiel gemaB den Fig. 5 und 6 ist eine Modification des in den Fig. 1 und 2 wiedergegegebe- 
nen Ausftihrungsbeispiels, bei dem zur Vereinfachung auf die Verarbeitung von StGtzwerten verzichtet worden 10 
ist Auf die Beschreibung der bereits in Fig. 1 und 2 angegebenen gleichartigen Baugruppen, die in Fig. 5 und Fig. 
6 mit gleichartigen Bezugszeichen dargestellt sind soli bei der Beschreibung dieses Ausfiihrungsbeispiels nur 
eingegangen werden, wenn es fur die Darstellung der Funktion erforderlich ist 

Durch die Modifikation des Filters gemaB Fig. 5 werden von dem Filter die Zieldaten eines Fahrzeugs in 
einem langeren Zeitintervall auch dann bestimmt, wenn sich das Ziel auf einer Zielbahn mit unterschiedlichen 15 
Geschwindigkeitskomponenten bewegt Diese Zielbahnabschnitte werden La. als "Legs" bezeichnet wobei 
angenommen wird, daB die zum jeweiligen Leg gehdrenden Geschwindigkeitskomponenten konstant sind 

Zur Detektion der Zeitpunkte, an denen das Zielfahrzeug ein Mandver, d h. eine Anderung seiner Geschwin- 
digkeitskomponenten vorgenommen hat, ist ein Mandverdetektor 76 vorgesehen, der mit der Peileinrichtung 87 
und der Steuervorrichtung 77 eingangsseitig verbunden ist Die Zeiten far detektierte Zieimandver werden in 20 
einem nachgeschalteten Mandverzeitspeicher 71 gespeichert, der einen ZShler 713 fQr die Anzahl detektierter 
Manoverzeiten und die Mandverspeicherplatze 71.0, 71.1 aufweist, wobei im Mandverspeicherplatz 71.0 die 
Anfangszeit T A o aller MeBzeiten 7i& dh. La. 75*1—0, gespeichert ist Dieser Mandverzeitspeicher 71 und der 
MeBzeitenspeicher 70 sind mit einer Zeitvergleichsschaltung 79 verbunden, deren Zeitausgfcnge 79.0 und 79.1 an 
die Koeffizientenschaltung 20 und an die Vergleichsschaitung 90 angeschlossen sind 25 

Im Zielvektorspeicher 10 sind fur die Geschwindigkeitskomponenten hach dem ersten Mandver und V y \ 
weitere Speicherelemente 16 und 17 vorgesehen, die ausgangsseitig einerseits mit der Koeffizienttenschaltung 
20 und andererseits mit dem Ausgabe-Ionterface 96 verbunden sind und die eingangsseitig ebenso wie die 
Speicherelemente 11 bis 14 durch die nunmehr erweiterte Torschaltung 95 und zur Initialisierung durch den 
Basisdatenspeicher 15 angesteuert werden. An die erweiterte Anzahl zu bestimmender Vektorkomponenten £0, 30 
?b, tfo Vxu %\ sind auBerdem die Rechenvorrichtung 80, die Akkumulationsschaltung 88, die Vergleichs- 
schaitung 90 und die Additionsstufe 94 angepaBt, was unmittelbar der erweiterten Anzahl von Verbindungslei- 
tungen zu entnehmen ist 

Die Koeffizientenschaltung 20 gemaB Fig. 6 weist weitere Eing&nge 21.4 und 22.4 fur die Geschwindigkeits- 
komponenten V x 1 und % 1 sowie einen weiteren MeBzeiteneingang 23.1 auf. An den Mefizeiteneingangen 23 und 35 
23.1 stehen als MeBzeiten die Zeitsignalwerte Tmm und Tm\ an. In Parallelschaltungen zu den Multiplizierern 
25.1 und 25.2 sind weitere Multiplizierer 253 und 25.4 vorgesehen, wobei der Muhiplizierer 253 eingangsseitig 
vom Eingang 21.4 mit der Geschwindigkeitskomponente 1 und vom Eingang 23.1 mit MeBzeiten Tm\ und der 
Multiplizierer 25.4 eingangsseitig vom Eingang 22.4 mit der Geschwmdigkeitskomponente %\ und vom Ein- 
gang 23.1 mit MeBzeiten Tm 1 angesteuert werden. Ausgangsseidg sind die Multiplizierer 253 bzw. 25.4 jeweils 40 
mit den Addierem 26.1 bzw. 26.2 verbunden. ^ 

Die Koeffizientenschaltung 20 enthait ferner weitere Multiplizierstufen 393 bzw. 39.6, die mit den Koeffizien- 
tenausgangen 56.1 bzw. 56.2 und jeweils mit dem Eingang 23.1 fur MeBzeiten Tmh verbunden sind Mit diesen 
Multiplizierstufen 393 bzw. 39.6 werden die zusatzlichen MeBkoeffizienten A/5 bzw. A/6 gebildet, die an den 
MeBkoeffizien tenausgangen 363 bzw. 36.6, die mit der Rechenvorrichtung 80 in Fig. 5 verbunden sind anstehen. 45 

In Fig. 5 und 6 ist zwar ein AusfQhrungsbeispiel fQr die Bestimmung der FilterkoeffizientenJ} bei einem 
einzigen Zieimandver dargestellt, jedoch durch gleichartige parallele Erweiterungen ist ein Filter zu realisieren, 
mit dem die Zieldaten eines Ziels bei einer beliebigen Anzahl j von Mandvern bestimmt werden. Darauf wird 
auch in mehreren der nachfolgenden Abschnitte hingewiesen. 

Mit dem modifizierten Filter gemaB Fig. 5 werden die Zieldaten eines Ziels bestimmt, das sich auf einer 50 
Zielbahn gemaB Fig. 7 bewegt Die Darstellung der Zielbahn in Fig. 7 ist auf ein kartesisches X K-Koordinaten- 
system bezogen, in dessen Ursprung sich die MeBposition X& Y& zu einem Zeitpunkt, L a. dem Anfangszeit- 
punkt der Messungen befindet Die Zielbahn beginnt wie dargestellt, zur Zeit Ta 0 an der Zieiposition Ai, Yo und 
das Ziel hat dabei die Geschwindigkeitskomponenten V* fy> Zum Zeitpunkt Ta 1 hat der Mandverdetektor 76, 
der das ZielgerSusch st&ndig einer Signal- oder Frequenzanalyse unterzieht, aus einer Signal- oder Frequenz&n- 55 
derung ein Mandver detektiert Das Ziel hat von diesem Zeitpunkt T A 1 an die Geschwindigkeitskomponenten 
9 kU V y \. Der Abschnitt def Zielbahn im Zeitintervall T A q, ta 1 wind auch als Leg 0 und der zum Zeitpunkt T A 1 
beginnende Abschnitt als Leg 1 beieichnet Bei der Detektion weitere Mandver wird die Zielbahn entsprechend 
ergSnzt, wie in Fig. 7 durch die gestrichelte Darstellung der Zielbahn fOr Leg 2 vom Zeitpunt T A 2 an angegeben 
ist Die zugehorigen Geschwindigkeitskomponenten w&ren dann Vx2 und ^2. 60 

Zur Bestimmung der Vektorkomponenten des Ziels auf den verschiedenen Legs werden nach der Anfangszeit 
Tao die Mandverzeitpunkte T AJ in den Mandverzeitspeichern 71 eingeschrieben. Die Detektion von Mandvern 
erfolgt immer erst nach einer durch die Konstruktion des Man6verdetektors 76 bedingten Anzahl von Messun- 
gen, also mit einer konstruktionsbedingten, zeitlichen Verzdgerurig. Das hat zur Folge, daB auch die Auswertung 
der Messungen durch das Filter erst nach eben dieser konstruktionsbedingten Verz6gerung erfolgen kann. 65 

Die Geschwindigkeitskomponenten V y bzw. % u %\ geben den richtigen Zielzustand paarweise nur in 
dem dem Leg zugehorigen Zeitintervall zwischen zwei Mandvern wieder. Daher gibt die Zeit zu Beginn eines 
Legs die untere Intervallgrenze und am Ende eines Legs die obere Intervallgrenze an. Bei dem Beispiel in Fig. 7 

11 
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ist fOr Leg 0 der Zeitpunkt T A0 die untere und der Zeitpunkt Ta i die obere Intervallgrenze, wohtngegen fur 
Leg 1 die untere Intervallgrenze durch Tax bezeichnet ist Die obere Intervallgrenze T A % wurde erst durch die 
Detektion eines weiteren ManSvers ermittelt Die dafur vorgesehenen Speicherplatze im Manoverzeitspeicher 
71 - sie sind in Fig. 5 nicht weiter dargesteUt - werden wie auch die Speicherplatze 71 A und 71.1 unmittelbar 
vorBeginndesMeBzyklusmitdenManoverzeitenrA/- » belegt. # 

In der Zeitvergleichsschaltung 79 in Fig. 5 werden die MeBzeiten 7m im MeSzeitenspeicher 70 nut den 
Mandverzeiten 7Vim Manoverzeitspeicher 71 vergUchen und den Intervallen, dh. den Legs toZidtahn, 
zugeordnet Jedem Leg ist ein Ausgang der Zeitvergleichsschaltung 79 zugeordnet Das heiBt, am Zeitausgang 
79.0 stehen Zeitsignalwerte Vm fur Leg 0 und am Zeitausgang 79.1 Zeitsignalwerte Tun fur Leg 1 an. Diese 
Zeitsignalwerte haben an jedem Ausgang den Wert Null, solange die MeBzeiten Tmi kleiner oder gleich der 
detektierten Mandverzeit Taj sind, wobei die Anfangszeit T A o in diesem Sinne\als Mandverzeit aufgebracht ist, 
Oder sie haben den Dif ferenzsignalwert zwischen MeBzeit Tm und unterer Intervallgrenze Taj. wenn die MeBzeit 
T^grdBer als die untere Intervallgrenze Taj und kleiner als die obere Intervallgrenze Taj+\ ist, oder sie haben 
den Differenzsignalwert zwischen oberer Intervallgrenze Taj^i und unterer Intervallgrenze Taj. wenn die 
MeBzeit TwigraBer als die obere Intervallgrenze T^+i ist Die Funktion der Zeitvergleichsschaltung ist for erne 
beliebige Zahl jvon Mane vern fur jeden Zeitausgang 79J durch die GL 15 

Tmu - max (0, min (Tab— Taj, Taj* t - Taj)) < 15 ) 

20 zu beschreiben. _ . . , ^ , . ,._ 

In der Koeffizientenschaltung 20 werden die Zeitsignalwerte 7TM;am Zeitausgang 79.0 mit den Geschwindig- 
keitskomponenten P„ V, und die Zeitsignalwerte Tm\ am Zeitausgang 79.1 mit den Gesc&windigkeitskompo- 
nenten fci, ?yi multipliziert, koordinatenweise zu den Komponenten Xo bzw. Y6 addiert, so daB sich daraus die 
erweiterten Zielpositionskoordinaten X b % ergeben. Unter BerOcksichtigung der zugehdngen gespeicherten 

25 MeBpositionskoordinaten Xm Ye ergeben sich dann analog den GL (3.1) und (3.2) die erweiterten Abstandskom- 
ponenten A* Ayizu 

A„ - to + % • Tm + Vyx ■ Tmx -Y B . I 16 - 2 ) 



10 



15 



30 



Mit diesen Abstandskomponenten A* A* werden die MeBkoeffizienten hu, . . .. A/4 gemaB GL (5) bestimmt, 
wobei sich die MeBkoeffizienten h l3 und A, 4 aus der Multiplikation der Me3koeffizienten An und h i2 mit den 
Zeitsignalwerten 7**1? an Eingang 23 ergeben. Die fflr die Ermittlung der Vektorkomponenten des urn die 
Geschwindigkeitskomponenten V xU t,\ erweiterten Zielzustandsvektor/) 

erforderlichen MeBkoeffizienten his und h t6 werden zusatzlich durch Multiplikation von An und A/a.mit den 
MeBzeiten Tmi an Eingang 23.1 berechnet und Ober die Ausgange 363 und 36.6 an die Rechenvomchtung 80 
Qbcrtr&crciL 

Wie sich aus den Bedingungen gemaB GL (15) fttr die Zeitsignalwerte Tm) an den Zeitausgangen ergibt ist 
derZeitsignalwert rammer dann idenrisch Null, wenn die MeBzeit Tm kleiner als die Mandverzeit 7^/ist.aber 
auch dann, wenn noch kein Man6ver detektiert worden ist Infolgedessen ist die Berechnung der erweiterten 
Filterkoeffizientenmatrix Fdurch die Rechenvomchtung 80 erst dann notwendig, wenn die Filterkoeffizienten- 
matrix Fttr eine aktuelle MeBzeit Tm bestimmt wird, die grdBer als die Mandverzeiten Taj ist ^Dazu wird die 
Anzahl j von Mandverzeiten Taj, die jeweils bis zur aktuellen MeBzeit Tm detektiert worden sind, in dem zur 
Adressierung der Mandverspeicherpiatze 71.0, 71.1 vorgesehenen Zahler 713 des Mandverzeitspeichers 71 
gezahlt und Qber den Steuerbus 78 an die Rechenvomchtung 80 Obertragen. Dadurch wird die zur Bestimmung 
von Koeffizienten h und /fOr eine Maximalzahl / von Zielmandvern ausgelegte Rechenvomchtung 80 derart 
gesteuert, daB jeweils der Grad der Matrizen tf und Fproportional der Anzahl /detektierter Mandver begrenzt 
ist Bei dem angegebenen AusfOhrungsbeispiel sind die Koeffizientenfolgen mandverbedingt, also jeweils urn die 
MeBkoeffizienten A^y+a und bzw. die Filterkoef fizienten / w+ 3 und erweitert 

In der Rechenvomchtung 80 werden gemaB GL (8) aus den MeBkoeffizienten hit bis A l6 die Filterkoeffizienten 
fn bis f i6 ermittelt, mit denen die Akkumulationsschaltung 88 angesteuert wird. Die Akkumulationsschaltung 88 
ist daher entsprechend der aufgrund einer maximal vorgegebenen Anzahl zu detektierender Mandver zur 
komponentenweisen Akkumulation der zu ermitteltennden Vektorkomponenten vorgesehen, dh. bei emer 
maximalen Zahl von /Mandvem sind 2/+ 4 paralle Akkumulationsstufen vorhanden. Die PrQfung und Verarbei- 
tung des erweiterten Ausgangsvektors Ap erfolgt grundsatzlich in einer Art und Weise, wie sei bereits bei dem 
Filter ohne zusatzliche Mandverdetektion gemaB Fig. 1 angegeben ist Lediglich in der Vergleichsschaltung 90 
ist far den Vergleich mit der Schwelle eine Vektornorm \\Ap || des Ausgangsvektors Ap 

zu ermitteln, die Komponenten des Ausgangsvektors Ap und nunmehr die Zeitsignalwerte Tmj an den Zeitsi- 
gnalausgangen 79.ybertlcksichtigt 
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Auch das Ausftthrungsbeispiel gemaB Fig. 5 ist durch die Eingabe von StOtzkoeffizienten h s zu verbessern. Die 
Stiitzkoeffizienten As ergeben dabei eine Verbesserung der Vektorkomponenten eines oder mehrerer Legs und 
sind entsprechend diesen Legs in die MeBkoeffizientenmatrix H eingeordnet Zu denjenigen Legs, zu denen 
keine weiteren stOtzenden Beobachtungen vorliegen, sind dann die zugeharigen Elemente der MeBkc>effizien- 
tenmatrix Hmit Nullwerten aufzufQllen. 5 
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